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It is shown that the Maxwell's equations for surface electromagnetic TM waves, propagating 
along the plane boundary between two nonlinear dielectrics with arbitrary diagonal tensor of 
dielectric permittivity, depending on |E|, can be integrated in quadratures. 



Surface electromagnetic waves, propagating along the boundary between two media with 
different optical properties were widely studied [1-7]. Nonetheless the integration of the 
corresponding equations presents difficulties because of their complicated nonlinear form. 
We will consider surface TM waves for the case of two dielectric media with the diagonal 
dielectric permittivity tensor e ik = diag^, e y , e z ), therewith e x , e y , e z are considered as 
arbitrary real functions of electric field strength E. 

Suppose the plane x = to be that of separation and plane electromagnetic wave with 
frequency u and wave vector k to propagate along z-axis. We will consider TM wave, 
assuming 

H= (0, 0); E= (E x ,0,E z ). 

Putting d t = —vd z , where v = ui/k = (5 c is the wave velocity, we write the system of 
Maxwell's equations in each medium: 

d x {e x E x ) + d z {e z E z ) = 0, (1) 

d x (H - f3s x E x ) = 0, (2) 
d z E x - d x E z = (5d z H. (3) 
Inasmuch at x = ± oo the field is supposed to vanish, from (2) we will find 

H = f3e x E x . (4) 

Assuming that d z = ik, E x = iA, E z = B, where A and B are real functions, from 
(1), (3) and (4) we deduce 

(s x A)' + ke z B = 0, (5) 

B' = kA(e x p 2 - 1), (6) 

where prime denotes differentiation with respect to x. Taking the functions e x (A 2 + B 2 ) 
and e z (A 2 + B 2 ) as given ones, we will admit that |E| 2 and e z can be expressed as some 
functions of e x , i.e. 

|E| 2 = A 2 + B 2 = I(e x ), e z (\E\ 2 ) = K(e x ) . (7) 
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Denoting e x A — /, from (5) and (7) we will find 



P f' 2 
1{£x) = 4 + k*K*(e x y 



where, according to (5) 



B = 



r 



kK{e x ) ■ 

Differentiating the equation (8) in view of (6), one gets 



dl , 



3-4 = (/ 2 )' + 



de 



1 



P 2 



K 



2f 2 _, 



3 4 • 



(8) 



(9) 



(10) 



The equation (10) admits the integrating factor Y(e x ) (see Appendix 1): 

de x 



\nY(e 3 



where 



XKAe r p * 



~k + 77k 



e*KX' 



(11) 



(12) 



Given £ x (0), the value of permittivity e x at E = 0, i.e. at x — ±oo, from (10) we deduce 
the quadrature 

Sx 

J Y(e x )dl(e x ) = f 2 X(e x )Y(e x ). (13) 

e*(0) 

The equation (13), thus, defines the function / = F{e x ), and finally one gets A = e~ l F(e x ). 

Solving the equation (8) with respect to /' and taking into account (13), we find the 
equation 



f = ^-e' x = ±kK(e x ) 



F 2 (e )iV2 



which is integrated in quadratures: 

kx = ± J 



dF{e x ) 



K{e x ) 



Ke x ) - ^ 



1/2 • 



(14) 



Here the sign ± is to be chosen in accordance with the domain x > or a; < 0. The integra- 
tion constants in (14) are to be found from the boundary conditions, which are equivalent to 
the continuity conditions for the functions / and /' at x = 0. The solution of the boundary 
conditions, mentioned, is defined by the form of the functions £a;(|E| 2 ) and e z (|E| 2 ). 
Let us consider a simple example assuming that 



£x — &x + &a;|E| 2 , By — Cl y + & y |E| 2 , E z — d z + b z \E\ 2 , 



(15) 



with a x , a y , a z , b x , b y , b z being constants. Then from (7) one gets the functions: 

I(e x ) = — , K(e x ) = q + — e x , 

o x b x 

with q = a z — (b z /b x )a x . Putting them into (11) we get (see Appendix 2) 



Y{e x ) = (q + ^e x ) 



2(3 2 e x + d 



d/b x ^Td 



(16) 



(17) 



.2f3 2 e x + d + y/D. 

where d = b z /b x — 1 and D — (1 + b z /b x ) 2 + Af3 2 q. Here we assume that D > 0. For 
simplicity we will now consider the case when the dielectric is isotropic, i.e. 



a + b\E\ 



(18) 



In this case one gets the integrating factor Y(e) = e. Putting it into (13) and taking into 
account that at x — > oo, / — > and e — > a, we find 



/ = m 

The equation (14) in this case reads: 

kx = ± J 



e(e 2 - a 2 ) 



\ 26(2 - (3 2 e)' 

3e 2 - (3 2 e 3 - a 2 
e(e - o)(2 - P 2 e)JP(7) 



de, 



(19) 



(20) 



where P(e) = (e + a) (3 e — 2 j3 e — a). One can rewrite this equality in the form: 



k x 



± 



a 

+ 2 



/-£- + ./■ 

J . m A J 



de 



P(s) ■' (e-a)jP(c) 
de . . af3 2 



+ 



r de a/r\ r 



de 



(2 - (Pe) y/Ptf 



Using the substitution 



cos2<y9, jx = V9 — 8 v, v = a{3 2 , 



A(3 2 A(3 2 

one can rewrite the last equality by means of elliptical integrals (see Appendix 3): 



kx = T[C F(ip, s) + J2 Cj n j(v?, n h s % 



(21) 



(22) 



(23) 



where 



Co 
Ci 

C 3 



2V2 



g2 = 2/ , 



e 2 ' 



8^ 



e(3-4i/+^)' 

£(3 + /i)' 
8y^2(l + i//2) 
^ (5 + A*) ' 



ni = 



2/i 



n 2 



™3 



3-4z/ + //' 
2/i 



3 + /i' 

2/i 
~5 + //' 



with f = y/'d + 4v + \x. 

The boundary conditions at x = 

e 1 E xl = e 2 E x2 , (24) 

and 

£*i = ^ 2 2 (25) 
for the case considered can be represented as follows: 

b 2 ei (ei - ai) (ei + ai) £ 2 {e 2 - a 2 ) (e 2 + a 2 ) 



6i /3i 2 £i-2 & 2 £ 2 -2 

(gi - gi + ri)(ei - qi - ti) (gg - q 2 + r 2 )(e 2 - q 2 - t 2 ) 
el (ei + ai) e| (e 2 + a 2 ) 



(26) 
(27) 

t 2 ^fc 2 -r u 2 ; 

where q t = 3/4(3? and n = ^9 - 8^/3?/ 4/3?, i = 1,2. 

The explicit form of relations (26) and (27) makes evident the existence of such param- 
eters &2/&1 or P for which the equations (26) and (27) are satisfied. 

Appendix 1: 

Multiplying the equation (10) by Y one gets: 

ds, • U ' \sl K + e x K) e\ ' 

The righthand side of the last equation can be written as: 



[f 2 XY]' - ^e'Y - fXY' - f 2 X'Y. 



'X 



Here X = l/e x — (3 2 / K + l/e x K. Equating the last three terms of the R.H.S. to zero we 
come to the equation 

r k> , 1 _ + «- 



Y XK 2K e x ' ' e 2 KX' 
which leads to the equation (11). 

Appendix 2: 
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Putting K(e x ) and X into (11) one gets 



\nY{e x ) = (b z /b x ) J^-dJ ■ 



K + IV 

with Q, — e x (1 — /3 2 e x ). Assuming that £> = (! + b z /b x ) 2 + 4/3 2 g > after integrating 



one finds: 



lny(e x ) = In | q + —e x | + [d/b x VD] In | — — — — = | . 

b x 2(3 2 e x + d + VD 



Appendix 3: 



For e = jjp + cos2(/9, one gets de = — svn2(p dip, and 



A(3 2 



In view of these expressions we get 

= 'Co [ , dV = -CoF{ V , s). 

J sJPie) J yjl - s 2 sinV 

Analogously one can find the other terms. 
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